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Abstract 

We investigate the possibility that the observed rotation of galaxies can be ac- 
counted for by invoking a massive baryonic disc with no need for non-baryonic 
dark matter or a massive halo. There are five primary reasons for suggesting 
this possibility. Firstly, that there are well known disc surface mass density 
distributions that naturally produce the observed rotation curves of galaxies. 
Secondly, that there are a number of rotation curve 'puzzles' that cannot be ex- 
plained by a massive dark matter halo i.e. the success of maximum disc fitting, 
HI gas scaling to the observed rotation, the disc/halo conspiracy and the inter- 
pretation of the TuUy-Fisher relation. Thirdly, recent 21cm observations show 
an almost constant HI surface density and a distinct 'cut-off' or edge to galactic 
discs. We explain this constant surface density in terms of either an optical 
depth effect or the onset of molecular gas formation and hence the possibility 
of considerably more gas existing in galaxies than has previously been thought. 
We suggest that the HI cut-off does indeed mark the edge of the galactic disc. 
Fourthly, there have recently been an increasing number of observations, most 
importantly 7-ray observations, that imply that Xco may be ten times higher 
in the outer Galaxy. This 'dark' gas may provide adequate mass to account for 
galaxy rotation. Finally, wc show that the additional baryonic mass required to 
account for the rotation of galaxies is just that required to reconcile observed 
baryons with those predicted by big bang nucleosynthesis without having to 
invoke a massive warm inter-galactic medium. We reconsider Mestcl's ideas 
about the collapse of isothermal and constant density spherical clouds and show 
that these can be simply used to siicccssfully model the rotation of galaxies. 
Mestel discs can also be used to straight forwardly explain the scaling laws of 
galaxies, particularly the observed relation between rotation velocity and radius 
and the oft used Tully-Fishcr relation. Thus the observed gross properties of 
disc galaxies can be explained by the monolithic collapse of baryonic gas into a 
rotationally supported disc. We discuss observations of the baryonic content of 
galactic discs and where sufficient 'hidden' baryons might be found to account 
for the rotation. 

1 Introduction 

Over the last 30 years dark matter has played an increasingly important role 
in our interpretation of Galactic, extra-galactic and cosmological observations. 
Circumstantial evidence for dark matter has gradually been accumulating from 



a wide variety of disparate observations; from locally observed large scale struc- 
ture (Doroshkevich et al., 2004) to the distant signals imprinted in the cosmic 
microwave background (Spergel et al., 2007). In this paper we want to take what 
some might describe as a 'blinkered' but, we hope, insightful look at one specific 
aspect of the dark matter problem. We discuss the requirement for a dark mat- 
ter halo around galaxies and specifically whether observations of galaxies could 
be explained more satisfactorily by a massive baryonic disc without resorting 
to non-baryonic dark matter at all. We do not address other evidence for dark 
matter. Our intention is to tackle a dilemma. Experiments have failed to detect 
dark matter particles, particle physics has failed to provide convincing dark par- 
ticle candidates and modifications to gravity do not sit comfortably within well 
tested general relativity. Perhaps one aspect of the dark matter problem, galac- 
tic rotation, can be explained by tangible well observed and familiar baryons 
without resulting to some new unknown form of matter. 

The almost constant velocity of rotation of galaxies at large distances from 
their centres (flat rotation curves) is most often quoted as the primary evidence 
for dark matter. In fact flat rotation curves in themselves provide no evidence 
at all because there is a disc surface mass distribution (see below) that naturally 
produces a flat rotation curve. What is important is not that the curve is flat, 
but the value of the asymptotic velocity and the large distances from the galactic 
centre that this velocity is measured. For a typical spiral galaxy, the asymptotic 
velocity at large distance from the centre is too great to be accounted for by the 
sum of the known baryonic components (stars, gas, dust) - hence the inference 
of dark matter. 

Before dark matter came into vogue astronomers were content to try and 
model galactic rotation using only the observed baryonic material. A seminal 
paper on this subject was written by Mestel (1963) (see also Hunter et al., 1984). 
The paper is quite extraordinary and brings into question the normally accepted 
and often quoted history of galaxy rotation, which firmly places the discovery 
of flat rotation curves in the 1970s, to early 1980s (Rogstad and Shostak, 1972, 
Roberts and Rots, 1973, Rubin et al, 1978) Q Mestel was trying to derive the 
surface mass density of two types of galactic disc that were distinguished by 
their observed rotation curves. The first had constant angular velocity (solid 
body rotation) and the second constant linear velocity. These two forms are just 
what are observed - many but not all galaxies have fiat rotation curves (Rubin 
et al, 1978). 

Mestel showed that these two rotation laws can arise from two quite simple 
disc surface mass density distributions. Constant angular velocity: 



^Though note Babcock's comment with regard to M31 that '...a very great proportion of 
the mass of the nebula must lie in the outer regions.', Babcock (1937). 
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and constant linear velocity: 

S(r) = ^ (2) 

r 

where E is the mass surface density, r the radial distance from the centre of 
the disc and Rq the size of the disc (r < Rq). (2) above is often referred to 
in the literature as Mcstcl's disc and is not exact (see below) because of the 
problem with the discontinuity at Rq. (1) above is what you would expect in 
the limit as a uniform density spheroid shrinks one of its axes to zero. This is 
important because Mcstcl's intention was not only to derive the surface mass 
distributions, but to also see if they were the viable result of the gravitational 
collapse of a cloud of gas - the then preferred mechanism for galaxy formation. 
He demonstrated that both arc acceptable solutions to the collapse problem 
provided the angular momentum distribution function (M(/i)) is preserved dur- 
ing the collapse i.e. that viscous transport is negligible - if so then (1) is the 
result of a collapsing uniform density cloud and (2) arises from a cloud in which 
the density 'decreases slightly' with increasing radius. 

M{h) is the function that specifies the amount of mass with angular momen- 
tum less than h. This function provides a test of Mcstcl's idea and ultimately 
an insight into the galaxy formation process. Crampin and Hoyle (1964) were 
the first to test Mestel's hypothesis that M{h) was the same for disc galaxies 
as it was for a uniformly rotating spheroid (the pre-galactic gas cloud). Using 
eight galaxies they confirmed Mestel's result. Essentially they measured the 
rotation and surface density as a function of radius to derive M{h). Freeman 
(1970) extended this work to a larger sample and came to the same conclusion 
as Crampin and Hoyle. That there was agreement between observation and 
Mestel's theory was in some ways fortuitous and in some ways insightful. It 
was fortuitous because at that time essentially only the inner solid-body part 
of the rotation curve was being measured. Also Freeman showed that M(h) 
for an exponential surface density distribution was very similar to the M{h) of 
Mcstcl's disc, thus importantly demonstrating that M{h) alone is not sufficcnt 
to determine the surface mass distribution of a self-gravitating disc. It was in- 
sightful because it provided something we no longer have - consistency between 
the observed distribution of matter (baryonic), how that matter rotates and a 
theory of how galaxies form. Today observations of disc galaxies have moved 
on, but maybe not forward. With a relatively high velocity flat rotation curve 
and a relatively low mass exponentially distributed baryonic disc M{h), at first 
sight, no longer provides a link to a plausible pre-galactic gas cloud. 

Our current theoretical ideas about galaxy formation are now quite differ- 
ent to those of Mestel. The hierarchical theory (White and Rees, 1978) places 
much more emphasis on galaxy merging as the process of large galaxy formation, 
rather than a single monolithic collapse. However, although there is clear evi- 
dence that galaxy merging does occur there is no clear evidence, for a galaxy like 
the Milky Way, that it is the major physical process at the root of its formation. 
Eggen, Lynden-Bell and Sandage (1962), using the metalicities and velocities 
of stars, were the first to present us with good evidence that the Galaxy is the 
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result of the gravitational collapse of a proto-galaxy - a slowly rotating cloud 
reduced to a thin disc within about 10^ years. A picture of early star formation 
in a spheroid (Population II) and subsequent star formation in a metal enriched 
thin disc (Population I) is still a viable explanation for at least some part of 
what we call the Milky Way galaxy. Based on the thickness of the Galactic disc, 
Toth and Ostriker (1992) have argued that less than 4% of the mass within the 
solar circle can have come from mergers and that this is inconsistent with the 
hierarchical galaxy formation models. This problem with the thinness of galac- 
tic discs has continued to provide a challenge to those who want to make galaxy 
mergers the major galaxy assembly mechanism (Navarro et al., 1994, Benson et 
al., 2004). 

If, as we believe, gravitational collapse still provides a viable explanation of 

the way a galaxy initially forms then might a Mcstcl disc be a viable explanation 
of the way it rotates, without the need for dark matter or alterations in the law 
of gravity? In this paper we want to reconsider Mestel's ideas. The posibility 
of a massive baryonic A^Icstcl disc is the first motivation for this paper. 

Our second motivation arises from some well known and long standing puz- 
zles regarding the rotation of galaxies. A rotation curve is normally modelled 
as the sum of the velocities of the presumed component parts of a galaxy: 

'^C = Stars + "^Gas + '^DM 

Each component's contribution to the circular velocity Vc is derived from the 
surface density of stars and gas and then the contribution of the dark matter 
is added so that the sum corresponds to what is observed. The dark matter is 
normally assumed to be distributed in an approximately spherical halo while 
the stars and gas arc in a disc. We, along with others before us (Pfenniger, 
Combes and Martinet, 1994 and references therein), identify four puzzles. 

1. Gas disc scaling - It was Bosma (1978, 1981) who first noted that to ac- 
count for rotation the required total surface density of matter in the spiral 
discs he studied was, in the outer regions, approximately proportional to 
the surface density of atomic hydrogen. In fact 'wiggles' in the rotation 
curves of some galaxies are directly reproduced by the inferred rotation 
predicted by the atomic gas alone (Fig. 1). There seems to be an intimate 
relation between the atomic gas and the dark matter. 

2. Maximum disc fitting - a problem with predicting the rotation due to 
the stars alone is the value of the mass-to-light ratio to use. Different age 
stellar populations and an imprecise knowledge of the initial stellar mass 
function can lead to quite different mass-to-light ratios. A solution has 
been to use the maximum disc model where the mass-to-light ratio of the 
stars is scaled so that the predicted rotation due to the stars fits the ob- 
served rotation curve in the inner regions. The deficit in rotation at larger 
radii is made up by the gas and dark matter (Palunas and Williams, 2000). 
The required maximum disc scaling is often consistent with expected stel- 
lar mass-to-light ratios and also with the stability of the stellar disc and 
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Figure 1: The lower dot-dashed hne is the rotation predicted for the gas alone. 
The upper black data points are the measured rotation curve. Note how the 
shape and even the 'wiggle' in the gas are reproduced in the observed rotation 
curve. The data shown are for NGC6946 - taken from Carignan et al. (2008). 



the formation of spiral arms (Fuchs 2002). Freeman (1992) showed that 
maximum discs provide an excellent fit to rotation curve data - there does 
not appear to be any need for dark matter in the central regions of a 
galaxy. 

3. The disc halo conspiracy - Given that over the inner region of a galaxy 
the stars and gas in a disc contribute significantly to the rotation while 
it is the spherical dark halo that does this in the outer regions, it is very 
strange that they 'conspire' to produce a flat rotation curve (van Albada 
and Sancisi, 1986). There is no known physical reason for why this might 
happen. 

4. The Tully-Fisher relation - The Tully-Fisher relation is one of the 
best-known scaling relations for galaxies (Tully and Fisher, 1977). It was 
originally proposed as a means of obtaining distances to galaxies for com- 
parison to their Hubble velocity distance, so that peculiar velocities could 
be derived. It has subsequently been used as a scaling relation in numeri- 
cal simulations to infer luminosity from measured dark halo rotation. The 
Tully-Fisher relation is between the observed velocity width derived from 
the 21cm line profile (Aw) and the luminosity (L). If galaxies are dark 
matter dominated then the line width is related to the mass (M) and ra- 
dial extent (R) of the dark matter while the luminosity is related to the 
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Figure 2: The HI surface density of THINGS galaxies scaled to the same maxi- 
mum radius - taken from Portas et al., (2008). 

baryonic material. The dark matter and luminous material again seem 
to have an intimate knowledge of each other. One interpretation of the 
Tully-Fisher relation is obtained by considering dynamical equilibrium: 



In this case E is now the surface brightness and Q the mass-to-light ratio. 
From the above the slope of the Tully-Fisher relation is predicted to be 
1/4 as long as Q and S are constant. The slope is observed to be, as 
predicted, approximately 1/4, but Zwaan et al. (1995) have shown that 
galaxies with a wide range of surface brightnesses obey the same Tully- 
Fisher relation - there seems to be some conspiracy such that as surface 
brightness decreases the mass-to-light ratio increases in exact proportion. 
There is also no known reason why this should be so. 

It is difficult to see how these four enigmas can be accommodated within 
our current theories of galaxies, galaxy formation and dark matter. 

Our third motivation for this paper comes from recent observations of atomic 
gas in galaxies. The THINGS team have been carrying out some of the most 
sensitive interferometric observations of galaxies using the 21cm line of atomic 
hydrogen (Walter et al., 2008). These detailed observations (Fig. 2) illustrate 
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some remarkable similarities for galaxies, which cover a large range in optical 
properties. The first intriguing thing is that they all have very similar HI col- 
umn densities of order lO^o-si cm~^, which remains reasonably constant with 
increasing radius over the entire disc (Bigicl et al. 2008, see also Minchin et 
al. 2004). The THINGS team's explanation of this is that at column densities 
above ~ 10^^ cm~^ the atomic hydrogen is converted into molecular hydrogen 
(H2). Interestingly, a column density of « 10^"'^ cm~^ is also just about the point 
at which atomic hydrogen becomes optically thick (see below) and so as well 
as molecules, there is also the possibility of more atomic gas hiding in the disc. 
The THINGS surface density profiles also show that there is an edge to the HI 
- a striking and rapid decline in the HI column density (Fig. 2). The THINGS 
team interpret this as the onset of ionisation of the HI by inter-galactic ultra- 
violet radiation. We will show below that the column densities of the discs of 
these galaxies are too high for such a sudden decline to be caused by ionisation. 
We alternately suggest that this drop in HI column density really does mark a 
physical edge of the galactic disc (Bland-Hawthorn et al. 1997). 

In summary, Mestel's ideas about gravitational collapse and the formation of 
a disc that has a flat rotation curve seem relatively simple, elegant and applicable 
to a galaxy like the Milky Way. There arc a number of disturbing puzzles 
about the way galaxies rotate that are difficult to explain if galaxy masses are 
dominated by a dark matter halo rather than just material in a disc. Finally, 
docs constant HI surface density provide a chic with regard to the distribution 
of baryons within the disc and does the sharp cut-off in the HI mark the edge of 
the disc? In this paper we explore the hypothesis that galaxies are massive gas 
discs, so massive that their rotation can be accounted for by baryonic material 
alone. 

Note that many of the issues highlighted in this paper have previously been 

discussed in Pfenniger, Combes and Martinet, 1994 (see also Pfenniger and 
Combes, 1994) and more recently by Hessman and Fiebart (2011). We hopefully 
provide further insight into essentially the same observational conundrums. 

2 Testing Mestel's ideas about gravitational col- 
lapse 

Mestel's calculations provide a simple, but eloquent analytical solution to the 
disc formation problem. We will now consider Mestel's conclusions in a little 
more detail and solve the equations numerically, that he looked to solve analyt- 
ically. First we will consider the assumptions that go into his model. 

1. A primeval gaseous sphere of 'roughly' uniform density (p) rotating with 
angular velocity contracts under its own gravity with each clement con- 
serving its angular momentum. So, the collapse occurs purely along the z 
direction i.e. parallel to the rotation axis. 

2. The turbulent and thermal motion of the gas make small contributions to 
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the momentum balance and so the centrifugal force acting on the gas is a 
measure of the radial gravitational field near to the galactic plane. 

3. The change in gravitational energy is dissipated in shocks and radiated 
away. 

4. No sub-condensations are formed and there is no further supply of kinetic 
energy by for example hot young stars. 

5. The collapse continues until a thin disc is formed and the infinitely thin 
disc approximation can be used. This requires that: 

which is true for a typical galaxy as long as the gas temperature is below 
about 10^ K (M ss 10^^ Mq, Rq « 30 kpc, c is the velocity dispersion of 
the gas). 

6. If the surface density changes with radius then so does the volume density 
and there will be a pressure gradient, but the force that this exerts is much 
smaller than that due to rotation and so it is ignored. 

7. The gravitational field perpendicular to the rotation axis changes so that 
centrifugal balance is maintained. 

Mcstcl's intention was to derive equilibrium distributions of mass and angular 
momentum of the disc given the conditions in the primeval sphere. The equi- 
librium rotation field is then: 



Q^ = -'^ (3) 
r ar 

The potential at r of a thin disc with size Rq and surface density S is given 

by: 



In each of the cases below we will assume a mass of M = 3 x 10^^ Mq and a 
radius of Rq = 30 kpc - a galaxy similar to the Milky Way. Let us first consider 
the collapse of the uniform density sphere. In this case the surface density is 
trivially: 
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Figure 3: The equilibrium rotation curves, calculated numerically, for the col- 
lapse of a uniform (blue) and power law density (n = 2) (black) spherical cloud. 
The blue dashed line is Mestel's (analytic) prediction of the constant angular 
velocity resulting from the collapse of the uniform density sphere. The 'Mestel 
disc' and exponential rotation curves are shown in green and red respectively. 
In each case the disc has a mass ofM = 3x 10^^ M0 and a radius of Rq = 30 
kpc, the exponential distribution has a scale length of 3 kpc. 




After considering the gravitational force exerted by an ellipsoidal mass distri- 
bution at a point on its equator as it becomes more flattened Mestel concluded 
that the resultant equilibrium disc would rotate with a constant angular velocity 
of: 

q2 _ SttGM 

In Fig. 3 we compare Mestel's prediction (blue dashed line) with that obtained 
by numerically integrating eq. 2 to get the potential and then using eq. 3 
to get the linear velocity fir (blue solid line). The discrepancy between the 
two is because Mestel's calculation is derived from the limit of a spheroid that 
shrinks one of its axis to zero and as such is not quite the same as doing the 
integration (numerically) over a finite disc. The conclusion though is essentially 
the same i.e. the collapse of a uniform density spherical cloud leads to a disc 
that rotates with approximately constant angular velocity. The surface density 
of the resultant disc is shown in Fig. 4, of particular interest (see below) is the 
almost constant surface density within the inner regions. 
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Figure 4: The disc surface density profiles, calculated numerically, for the col- 
lapse of a constant density (n = 0) (blue) and power law (n = 2) (black) 
spherical cloud. 'Mestel disc' and exponential surface densities are shown in 
green and red respectively. Each profile is for a disc mass oi M — 3 x 10^^ M0 
and radius of Rq — 30 kpc, the exponential distribution has a scale length of 3 
kpc. 

The other interesting case is that of 'Mestel's disc' and what density spheroid 
gives rise to it? In his paper Mestel comments that a 1/r surface density distri- 
bution results from the collapse of a spheroid that has a 'slowly varying density', 
he does not quantify this. We have experimented with initially spherical power 
law mass density distributions to get the best chi squared fit to a 1/r surface 
density law, there being no exact fit using a power law. To two decimal places 
we find that n = 1.76 gives the best fit to a 1/r surface density distribution 
given the spheroidal mass density distribution p{w) = po{w / Rq)^ , where w is 
the radial co-ordinate of the spheroid. This is very close to the n = 2 case which 
is the well known density distribution for the singular isothermal sphere i.e the 
collapse of a singular isothermal sphere gives rise to a surface density distribu- 
tion that is very close to a Mestel disc {^(r) = Eo(r/i?o)"° with Ug —1). In 
Fig. 4 we show the resultant surface density distribution (n — 2) along with 
that of the exact (1/r) Mestel disc (black and green lines), they are very similar. 

As for the constant density spheroid case we can now use equations 1 and 2 to 
obtain the predicted rotation curve. Mestel's conclusion with regard to constant 
linear velocity was derived by simplifying a series expansion and assuming that 
the disc size tends to infinity. In Fig. 3 we show the result of the numerical 
solution of the predicted rotation of a finite disc both for the case where the disc 
surface density varies as 1/r and that obtained for the collapse of a sphere with 
density power law index n = 2. Mestel's assertion of constant linear velocity is 
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almost, but not quite, correct for the finite disc as the velocity increases from 
about 200 to 230 km over the radius of the disc (Fig. 3 green line). The 
predicted Mestel constant disc velocity (i/GM/J?05 see below) is actually about 
211 km s^^. The predicted velocity for the collapsed n = 2 power law is similar 
to the 1/r case in that it produces almost constant velocity (Fig. 3 black line). 
Our conclusion is by no means new, but because of its simplicity it is worth 
repeating: we can achieve an approximately constant linear velocity disc from 
the collapse of a singular isothermal sphere. 

Before the advent of dark matter the dominant mass component of galaxies 
had been taken to be the stars. The challenge then was to model the rotation 
using essentially the stars alone. The stars were observed to follow an expo- 
nential surface brightness distribution (Freeman 1970) and so alone they could 
not account for the rotation. For completeness in Fig 3 and 4 we show the 
corresponding rotation curve and surface density for the exponential disc. 

The specific cases of constant angular velocity and constant linear velocity 
are of particular interest for those interested in the origin of galaxy rotation 
curves. Ignoring the fine structure, many (all?) rotation curves can be mod- 
elled as combinations of constant angular and linear velocity components. Dwarf 
galaxies often have continuously rising rotation curves (constant angular veloc- 
ity) while galaxies like the Milky Way have almost flat rotation curves (constant 
linear velocity) over most of the disc. Other galaxies have some combination of 
the two. Following Mcstcl's reasoning, the implication is that the original gas 
clouds that galaxies formed from, assuming that this was monolithic collapse 
have, to varying degrees, constant density cores surrounded by more diffuse 
isothermal spheres. 

Let us now consider some speciflc cases. We have taken data points from 
the rotation curves of three galaxies studied as part of the THINGS survey (de 
Blok et al., 2008). The galaxies were chosen to illustrate the resultant rotation 
of discs that collapse from singular isothermal spheres, constant density spheres 
and some combination of the two. NGC7331 can be modelled as a singular 
isothermal spherical cloud {p{w) cx w^^) of mass M = 2.8 x 10"'^'^ Mq and size 
Rq = 30 kpc. With this density law the cloud collapses to give a surface density 
distribution (S(r)) that is approximately that of the 'Mestel disc' (a r~^) and 
hence an almost flat rotation curve (Fig. 5). NGC2976 is a much less massive 
galaxy and has a rising rotation curve. This can be modelled as resulting from 
the collapse of a M = 4.5 x 10^ Mq constant density cloud of size i?o = 3 kpc 
(Fig. 5). NGC2403 requires a combination of a constant density core (size« 2.5 
kpc) and p{'w) cx outer regions. The total mass is M = 6.0 x 10^" Mq and 
size Ro = 20 kpc, 14% of the mass is contained within the constant density core 
(Fig. 5). 

The rotation curves of the above three galaxies can all be explained as result- 
ing from the collapse of gas clouds with the most simple and plausible density 
distributions. Small clouds are of almost constant density while large clouds 
become more diffuse in their outer regions. As far as we are aware all galaxies 
have constant angular velocity as you approach their centres, which must arise 
from an initially approximately constant density core in the primeval cloud. The 
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Figure 5: The equilibrium rotation curves fitted to observational data for three 
galaxies taken from the THINGS survey. The origin of NGC7331 is almost 
entirely due to the collapse of a singular isothermal sphere gas cloud, NGC2403 
requires a constant density core within the isothermal sphere while the origin 
of NGC2976 is a constant density cloud. 



collapse of this constant density core leads to an almost constant surface density 
(Fig. 4). This measurement of a constant density at the centres of galaxies (de 
Blok et al. 2003) has been a long standing thorn in the side of those who model 
galaxies using cold dark matter halos. Cold dark matter haloes are predicted 
to have steep 'cuspy' density profiles (Navarro et al. 1996). 

3 Scaling relations 

If disc galaxies are primarily Mestel discs we can very easily predict what the 
fundamental scaling relations between their physical characteristics should be. 
It is clear from the above that the three quantities that define a Mestel disc 
are the radius {Rq), total mass (M) within Rq and the velocity («„) at Rq. 

Where for stability u,„ = \J~^^i so just two of these three parameters are 
independent. We have used the extra-galactic database LEDA to extract data 
for a large (5174) sample of galaxies. Galaxies were selected such that they have 
angular diameters (r25) greater than 1 arc min, an I band magnitude, a 21cm 
line flux, a measure of the maximum rotation velocity of the gas and finally that 
they have a morphological type oit> 2 (Sab or later). 

Relatively large angular sized galaxies were chosen so that hopefully they 
will have been accurately measured, particularly their sizes. I band magnitudes 
and 21cm line fluxes were required to calculate stellar and gas masses. We can 
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Figure 6: The velocity/radius fundamental scaling relation using 5174 galaxies 
from the LEDA database. The solid blue line is the best fitting line of slope 
unity. The two dashed lines are for densities a factor of three lower (upper 
line) or higher (lower line) than the best fitting line. The black line is what is 



expected if galaxies are assembled directly as discs so that c>c i?. 
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equate the LEDA 'maximum rotation velocity of the gas' directly to Finally, 
by using t > 2 we hopefully restrict our analysis to galaxies dominated by their 
disc. 

It is clear that we have independent measures of the velocity and size and 
so it is the mass of the galaxy that becomes the dependent variable. The only 
remaining problem is how does the size given in LEDA (r25) equate to the size 
Rq used in our model (Note: Rq is the physical size of the disc), which we are 
assuming is the same as the cut-off observed at 21cm (see below). During their 
analysis of THINGS galaxies Portas et al. (2008) quote the measurements of 
Broeils and Rhee (1997), who give a value of rif j /r25— 1-7 (±0.5) measured at 
a column density of 1.25 x 10^" cm~^. This column density seems to agree very 
well with the cut-off seen in the surface density profiles of Fig. 2 and is the value 
we will use to scale the optical to the HI disc cut-off radius. Thus, we now have 
the two independent variables we require, Vm and Rq, they are plotted against 
each other in Fig. 6. 

Let us assume that isothermal spheres give rise to Mestel discs that produce 
the flat rotation curves we observe. If all galaxies initially form from similar 
density gas then there should be a mass size scaling relation M (x Rq. Using 
the Mestel disc equation for velocity and the substitution for M leads to the 
prediction Rq (x Vm- Although true for Mestel discs this Vm/Ro relation is not 
universally true for all disc surface densities. For example from Fig. 3 it is clear 
that the maximum rotation velocity of the exponential disc is much higher than 
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that predicted by a Mestel dise or by the collapse of the uniform sphere. So it 
is important to note that the above scaling relations only apply to discs that 
result from the collapse of spheres with < n < 2. On Fig. 6 we have drawn 
the best fitting line of slope unity and it is clear that it fits the data very well. 

However, the measured gradient of the best fitting line is 0.8 and not 1.0, but 
there is an obvious explanation of this discrepancy. Lines of slope unity in this 
plot correspond to lines of constant initial density in the proto-galactic spheroid 
(po) and so the fit is dependent on the relative numbers of galaxies of different 
densities along the line. The range of densities is actually remarkably small, 
the blue dashed lines on Fig. 6 illustrate a change in initial proto-galactic cloud 
density of only a factor of 3 from the mean, with most galaxies being constrained 
within these two dashed lines. 

If galaxies were assembled directly as discs then we might expect a different 
relation between mass and size i.e. M oc and consequently that Rq oc as- 
suming in this case that the surface density remains constant. This relationship 
is also plotted on Fig. 6 (black line) and seems to be at odds with the data. It 
appears that galaxies retain a memory of their spheroidal origin. If isothermal 
spheroids collapse to Mestel discs then we have M = 4:npoRo — SttSo-Rq ^ re- 
lationship (So = 2poRo) between the measured surface density of the disc and 
the density of the cloud it formed from. 

We can now define a 'typical' galaxy using the LEDA data. The mean 
maximum rotational velocity of the sample galaxies is v„i ~ 144 km s""'^, mean 
radius Rq ~ 20 kpc and resultant mean mass M/jot = 1.6 x 10^^ Mq. The 
'characteristic' density (pg) of the cloud that this 'typical' galaxy formed from 
is 7 X 10~^^ kg m^. If the current value of the baryonic mass density of the 
Universe is ~ 2 x 10"^^ kg m""^ (Spergel et al. 2007) then is achieved at a 
redshift of about 14 - not an unreasonable redshift for galaxy formation. 

We now need to consider how much mass is required to be consistent with the 
observed velocities and sizes of our sample galaxy. The predicted (dynamical) 

2 p 

mass of each galaxy is simply calculated using M = Substituting for 

Vm or Rq in this equation leads to Rq oc M^/"^ and Vm oc M^/^ (we also have 
So oc M^/^). We suggest that the latter relation {vm oc M^/^) is the origin 
of the Tully-Fisher relation (see below). In Fig. 7 we show these two scaling 
relations for our LEDA sample. Again variations in initial density will produce 
scatter in the plots, but they are both extremely well fitted by a straight line 
with exponents of 0.33 and 0.34, for the velocity/mass and radius/mass relations 
respectively. This is remarkably close to expected value of 1/3 for Mestel discs. 

We have also used the LEDA I band magnitudes and 21cm line fluxes (/) to 
calculate the stellar and atomic gas masses of our sample galaxies. Stellar masses 
have been simply calculated assuming Sun like stars and M® = 4.08. Atomic gas 
masses are calculated using the standard formula Mhj = 2.4 x lO^d^^p^ / fdv. 
In each case we have used the LEDA recession velocities and a value of Hq = 72 
km Mpc"-^ to get distances. With a mean stellar mass of 1.8 x 10''^'' Mq 
and atomic gas mass of 0.7 x 10^° Mq there is not a huge discrepancy between 
the baryonic mass we can account for and the mass required for rotation (< 
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Figure 7: The fundamental scaling relations of galaxies using 5174 galaxies from 
the LEDA database. Top - the mass size relation. Bottom - the mass maximum 
velocity relation. In each case the blue dashed line is a linear least squares fit 
to the data. 
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Figure 8: The blue histogram shows the distribution of values of the ratio of, 
dynamical mass to mass in stars and atomic gas. The red histogram shows the 
distribution of the value of k that is needed to multiply the atomic gas mass by 
so that the rotation can be accounted for. 



Mnot >— 1-6 X 10"'^^ Mq). The mean ratio of dynamical mass to stellar and 
atomic gas mass for the whole sample is w 7 (Fig. 8). This scaling factor of 
7 is consistent with the cosmological 'missing mass' problem which measures 
^Matter ~ 0.24 and ilsaryon ~ 0.04 leading to = 6 (Spergel et al., 

2007). 

As we will see in the following sections we will suggest that the mass discrep- 
ancy between what we can 'see' {Mstar + Mhi) and that measured via rotation 
{Mfiot)i is 'hidden' in a previously undetected gas component. To assess just 
how much gas is 'missing' we also calculate the factor k which multiplies the 
atomic gas mass to fully account for the dynamical mass i.e. k — ^^"°Mjff^'°'' ' 
its mean value for the LEDA sample is 27 and its distribution of values is also 
shown in Fig. 8. 

Given the gas disc scaling relation (Bosma 1978, 1981) and that the observed 
surface density of stars appears to be exponential we require an additional mass 
component, larger than that of the stars, but distributed like the gas to account 
for the rotation. As a demonstration that models like this are possible and 
stable Kalberla and Kerp (1998) demonstrate that a thick (« 4 kpc) disc of 
mass w 5 X 10^^ Mq can account for the rotation of the Milky Way (see also 
Kalberla, 2003, 2004, Kalberla et al. 2007 and Revaz et al. 2009). 

In the following sections we will consider observations of the mass of gas 
in galaxies, how uncertain these might be and whether this uncertainty can 
account for the value of k and ultimately the dynamical mass of galaxies. 



16 



4 Observation of gas in galaxies 



4.1 Atomic hydrogen 

As noted above the THINGS team measure an approximately constant atomic 
hydrogen surface density for galaxies (Bigiel ct al., 2010). The value of the 
observed HI column density (10^° < N < 10^^ cm~^) is intriguing. The optical 
depth of atomic hydrogen (Kerr, 1975) is: 

5.5 X IQ-^'^N 

For Ts = 100 K , cr = 10 km s~^ clouds with column densities (iV) of about 
10^^ cm~^ or above are opaque, r is inversely proportional to T and Ts =100K 
is an upper limit on the cool atomic gas temperature measured by Dickey et 
al. (2009). Almost all inferences made about HI in galaxies assume that it is 
optically thin. 

lonisation of the HI will occur (re-combinations=ionisations) when N k. 
where </> « 10 photons cm~ s~ is the intensity of the ionising background and 
a = 3 X 10^ ^^Tg''''''' cm'' is the hydrogen recombination coefficient (Maloney, 
1993). For a disc with a gas z scale height of (3 the above equation becomes: 

N ft! \J For a gas scale height of /3 = 500 pc and again a temperature of 

100 K, the column density at which ionisation becomes important is N k. 10^^ 
cm~^. At column densities lower than this there is an increasing ionisation 
fraction. The typical explanation of the sharp cut-off to the HI profiles shown 
in Fig. 2 is that this is due to the ionising background. But, as noted by Portas 
et al. (2008), the drop in column density starts at ~ 10^° cm~^, which is too 
high for the HI to be substantially affected by ionisation. Hence our suggestion 
that this really does mark the edge of the disc. 

The predicted range of observed column densities due to ionisation and op- 
tical depth (10^9 <N < 10^1 cm'^) and the THINGS result for other galaxies 
(10^° < N < 10^^ cm~^) can be compared with observations of our Galaxy. 
Stark et al. (1992) made a map of the HI column density across the sky and 
as expected the measured column densities are somewhat restricted (Fig. 9). 
There is quite a sharp cut-off at N = 10'^^ cm~^ and then it just slowly falls 
away to 10^^ cm~^. Most lines of sight lie in the predicted range 10^°~^^ cm~^. 
The Stark et al. survey column density limit is « 2 — 3 x 10^^ cm~^. Thus the 
cut-off at low column densities (Fig. 9) appears to be real rather than an arte- 
fact of the observational set up. One interpretation of the narrow range of 21cm 
column densities observed for both our Galaxy and THINGS galaxies is that 
the highest column density sight lines are optically thick and hiding substantial 
amounts of atomic gas. 

An alternative way of measuring the HI column density along a line of sight, 
that does not rely on 21cm emission, is to use stellar spectra. Abundant atomic 
hydrogen should cause strong absorption particularly in the Lyman q line. 
Bohlin et al (1978) have mapped the HI distribution using Lyman a absorption 
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Figure 9: The distribution of atomic hydrogen column densities from the Stark 
et ah (1992) observations of our galaxy (black). The red line shows the column 
densities derived from ultra-violet observations of nearby stars (scaled to give 
the same peak value). 



lines seen in the spectra of 100 stars observed by the Copernicus satellite. In- 
terpretation of these lines is not straight forward because the limiting distance 
for these ultra-violet observations is dependent on the dust extinction, which 
is itself dependent on the HI column density. So the densest gas regions are 
un-observable. Bohlin et al. placed a limit of E(B-V)< 0.5 for sight lines to 
the stars they observed, which translates, into a HI column density of less than 
2 X 10^^ cm~^. The distribution of column densities sampled by the stars is also 
shown in Fig 9., where it is compared to that derived from the 21cm emission 
data. It is clear that the distributions are very different making it difficult to 
say if either is a true reflection of the real atomic column densities through the 
disc. What is interesting is that the stellar data has increasing numbers of sight 
lines with increasing column density while the HI data has decreasing numbers 
with increasing column density. One can speculate whether this trend for the 
stellar data would continue, and hence make the inference of much more atomic 
gas, if it was not arbitrarily cut-off at just over 10^^ cm^^. Bolin et al. also, 
rather interestingly, point out that in a number of cases the 21cm and Lyman 
a measurements do not agree on the measured column density. If optically thin 
the 21 cm observations should sample a longer HI column than the distance to 
the star so -^Ly ^ -^21cm' ^^^^ ^'^^ always the case. 

How much atomic gas could there be in galactic discs without contradicting 
the Stark et al. and THINGS observations? Let's return to our 'typical' Mestel 
disc. In his paper Mestel (1963) also solves the problem of the predicted surface 



18 



density distribution of a truncated disc that has constant Unear velocity finding 
that: 



This now has zero column density at and beyond the truncation radius Rq. We 
will use this truncated disc to model the opacity effect and what it predicts. 

By requiring that v.m=211 km s^^ and that i?o=30 kpc for our 'typical' 
galaxy we have a total disc mass of M = 3 x 10^"'^ Mq, which wc now assume is 
all in the form of atomic hydrogen distributed to an exponential scale height of 
0.5 kpc (Bohlin et al. 1978). We have constructed a simple numerical radiative 
transfer model of a Mestel gas disc using the above parameters, including the 
effects of 21cm opacity (eq. 5 above) so that we can compare with observations. 
The predicted observed (face-on) radial 21cm column density distribution is 
shown in Fig 10 along with the actual column density. The predicted almost 
constant observed column density has about the correct value of 10^"'^ cm~^ 
(compare with Fig 2) while the real column density is more than an order 
of magnitude higher over most of the disc. This HI disc, massive enough to 
account for its rotation, has an almost constant observed surface density of 
hydrogen consistent with observation. There is however a 'catch'. To achieve 
the necessary opacity we had to use a gas temperature of lOK (and a =10 km 
s~^). Thus to be able to account for the rotation using atomic hydrogen alone 
we need it to be cold. Can cold atomic gas be present in galactic discs? There 
are three types of observation that provide us with information on the physical 
state of the atomic gas. 

The first and most familiar is the emission line surveys such as THINGS. The 
difficulty with these is that the emission is independent of the gas temperature 
when it is optically thin and proportional to the temperature when optically 
thick. Thus observations of the emission do not lead directly to a measurement 
of the column density unless the gas is assumed to be optically thin, which is 
what is almost universally done. Returning to the Stark et al. observations of 
the distribution of HI column densities on the sky (Fig. 9) these can actually 
be interpreted in two ways. Previously the HI has been assumed to be optically 
thin so the observed distribution of column densities would be regarded as a true 
measure of the column densities along each line of sight through the Galaxy. 
This is not so if these sight lines are optically thick. Using our above radiative 
transfer model of a typical galaxy we can predict what we would observe from 
about the position of the Sun in the mid-plane. For gas with Ta = 100 K 
km s^^ the observed column density would be about lO^'' ''' cm^^ in all direc- 
tions independent of the actual column density. For Ta = 1000 K km s~^ the 
observed column density increases to about 10^^'^ cm~^. Thus an alternative 
interpretation of the Stark et al 21cm data is that it is really a measure of Ta 
along the line of sight in the optically thick regime and large amounts of atomic 
gas could be hidden there. 

Secondly, 21cm absorption studies generally target background continuum 
radio sources so that absorption by hydrogen in our galaxy can be studied 
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Figure 10: The real (black) and observed (red) atomic hydrogen column density 
profile of an optically thick face-on Mestel disk. In this case the temperature of 
the gas is lOK, it has a velocity dispersion of 10 km s~^ and a z scale height of 
0.5 kpc. Note for comparison with Fig. 4 that 10^^ cm~^ « 8 x 10^ kpc~^. 



(Dickey et al. 2009, Heiles and Troland, 2003, Kolpack et al. 2002, Garwood 
and Dicky, 1989, Dicky et al. 1978). These observations are primarily con- 
cerned with absorption 'features' and as such they rely on the premise that the 
hydrogen is contained within discrete clouds - it is clearly much more difficult to 
detect hydrogen extended both spatially and in velocity. The surveys do provide 
useful physical parameters for the gas in clouds, which is normally divided into 
cold and warm components. Typically the cold component is observed to have 
a spin (kinetic?) temperature of about 50K with some less than 25K, while 
the warm component has T > 500K. Absorption by the warm component is 
generally at larger distances above the Galactic plane than the cold component. 
Interpretation of the absorption lines is fraught with difficulties: lines may be 
blends from different clouds, the on source absorption spectrum may contain a 
contribution from HI emission in the beam, baselines are difficult to measure 
and to measure optical depths you need an emission measure that is determined 
off-source. Data is also often difficult to interpret because authors give proper- 
ties of the clouds, but not the integrated optical depth along the line of sight. 
It is the cold component (eq. 5) that gives rise to high optical depths, the warm 
component's optical depth being generally pretty low. 

An exensive study of 21cm continuum source absorption through the gas disc 
of our galaxy was carried out by Kolpak et al. (2002). Their stated sensitivity 
gives a minimum detectable column density of 7 x 10^^ cm^^ though this is 
in 'clouds' and so provides only an estimate of the minimum column density 
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(and opacity) along a line of sight. Kolpak et al.'s objective was to calculate 
a radial 21cm opacity (r) profile for the Galaxy. They reach four important 
conclusions. Firstly, that within 3-8 kpc from the Galactic centre (they have 
few sight lines within 3 kpc) the HI opacity lies between t =0.8 and 1.5 with 
some sight lines having opacities as high as 3.5. The inferred column densities 
of hydrogen depend on the temperature and velocity dispersion along the line of 
sight, but these observed opacities are not at all inconsistent with our optically 
thick 'typical' galaxy. Kolpak et al. say that there are saturated lines at all radii 
and presumably these are given a minimum r consistent with this saturation. 
Secondly, low opacities arc measiired beyond 8 kpc, but at these radii they say 
opacities are very uncertain because of the uncertainties in the rotation curve. 
Thirdly, there is a good correlation of high HI opacity with high molecular gas 
emission, as might be expected if cold HI eventually loads to the formation of 
molecular clouds. Fourthly, there is quite strong disagreement about the opacity 
distribution compared to a similar study carried out by Garwood and Dickey, 
(1989), which further illustrates the difficulty of interpreting absorption line 
studies. 

The third and final way of measuring the quantity of neutral hydrogen is 

by HI self-absorption. HI self-absorption occurs when there is gas with a low 
excitation temperature lying infront of a hotter cloud. This of course requires a 
rather special geometrical alignment, but in the plane of the galaxy this seems 
to happen quite often (Gibson et al., 2005). Strong HI self-absorption again 
requires significant density enhancements (clouds) and a large scale background 
at the right velocity and/or very cold gas. This self-absorption does not provide 
a strong constraint on the total amount of atomic hydrogen present, but Gibson 
et al., (2005) do say, after carrying out an extensive self- absorption study of the 
Galactic plane that '...faint HI self-absorption is detected virtually everywhere 
that the HI emission background is sufficiently bright'. 

There are of course quite different problems associated with discerning the 
mass of atomic hydrogen in our Galaxy compared to other galaxies because we 
reside amongst the hydrogen. Our current understanding of the Milky Way 
has benefited greatly from observations made of other galaxies, particularly our 
nearest neighbour M31. Recently Braun et al. (2009) have found what they 
describe as 'ubiquitous HI self-opacity features' over the face of the disc of M31. 
They say that in some cases there is an order of magnitude difference in the 
HI column density inferred using the normal optically thin assumption, but 
that currently they account for only about a 30% increase in the previously 
determined HI mass. This they say though must be a lower limit because of the 
geometry requirements of HI self-absorption stated above and the blending of 
lines. They find gas spin temperatures of 20-60 K for the opaque gas, consistent 
with those required to give high opacity in our 'typical' model galaxy. It is also 
worth noting that they find a strong truncation of the HI disc at 28 kpc (at 
5 X 1019 cm-2). 

Perhaps the most convincing observation consistent with the idea that there 
is both an extensive and massive cold HI component to galaxies comes from the 
observations made by Knee and Brunt (2001). We quote from the abstract of 
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their paper - 'Cold HI is difficult to observe, and therefore our knowledge of its 
abundance and distribution in the interstellar medium is poor. The few known 
clouds of cold HI are much smaller in size and mass than typical molecular 
clouds. Here we report the discovery that the HI supershell GSH139-03-69 is 
very cold (lOK). It is about 2 kpc in size and as massive as the largest molecular 
complexes. The existence of such an immense structure composed of cold atomic 
hydrogen in the interstellar medium runs counter to the prevailing view that 
cold gas resides almost exclusively in clouds dominated by molecular hydrogen'. 

4.2 Molecular hydrogen 

Another possible substantial component of the interstellar medium in the discs 
of galaxies is molecular gas, primarily in the form of molecular hydrogen. The 
direct observation of molecular hydrogen presents substantial problems because 
it has no permanent dipole, so there is no pure rotational spectrum, which for 
molecules would normally produce emission in the ground based mm window. 
H2 rovibrational lines can be detected at wavelengths of a few micron, but these 
arise from gas at a few lOOOK and so it does not correspond to the cold molecular 
material (10-40K) that is known to reside; in the discs of galaxicis. Generally the 
quantity of molecular hydrogen has been inferred from the relatively abundant 
CO molecule that does have a dipole moment and so emits strongly (Dickman 
1978). 

The exact relation between the abundance of H2 and CO has been a bone 
of contention over the years, with a current best buy value of Xqo ~ 2.0 x 10^*^ 
cm~^ K^^ km^^ s where Xco — N{H2) / Ico, N{H2) is the molecular hydrogen 
column density and Ico is the brightness temperature of the J — 1 — emission 
line at 2.6mm (Burgh et al. 2007). Historically confidence in both the constancy 
of Xco and its precise value comes from three independent measures that concur 
(Solomon and Barrett, 1991): 

1. The correlation of CO intensity with extinction (Ay) determined by star 
counts. The H2 column density is then inferred from the Nh/Av relation. 

2. By modelling the 7 ray flux from the interaction of cosmic rays with 
hydrogen molecules and comparing to the CO intensity. 

3. Virial analysis using cloud sizes and line widths to get the total molecular 
mass and then comparing this mass to the CO intensity. 

Contrary to the above there are two quite logical reasons why we might not 
expect a good correlation between CO and H2 and over more recent years these 
have been investigated. Firstly, the quantity of CO compared to hydrogen 
should depend on how metal enriched the gas is and so Xqo should depend 
on metalicity. Estimates of the dependence range from Xco oc (Israel 
2000) to Xco oc Z-^ (Wilson 1995, Boselli et al. 2002). Given the metalicity 
gradients found in galaxies (Magrini et al. 2011) it is easy to see how the former 
(oc Z~'^-^) metalicity dependence leads to the conclusion that Xco could be a 
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factor of 10 higher in the outskirts of galaxies. For example using a sample of 
29 SINGS galaxies Draine et al. (2007) have derived both HI and H2 masses. 
Using their values we obtain a mean value of < k >= {Mhi + ) I^hi = 3.2. 
This is based on a constant value of Xqo^ which if made 10 times larger would 
give values of < fc > just about that required to account for galactic rotation 

In support of the above are observations of CO absorption in cold (« lOK) 
clouds in our galaxy at a radial distance of about 12 kpc (Lequeux et al. 1993). 
At this distance Lequeux et al. say that the mass of H2 is already 4-5 times 
large than HI (fc = 5 — 6) and that it will only increase with increasing distance. 

Secondly, models indicate that the photodissociation of CO is more sensitive 
than H2 to the strength of the interstellar UV radiation field (van Dishoeck and 
Black 1988, Kopp et al. 2000). Thus in different galactic environments we might 
expect different ratios of CO to H2. Recently Ade et al. (2011) using PLANCK 
data have used just such an explanation to interpret an observed breakdown in 
the correlation between optical depth and hydrogen column density (between 
about 8 X 10^° and 3 x 10^^ cm~^). They suggest that the CO is being pho- 
todissociated and so it no longer becomes a good tracer of H2. They describe 
the unaccounted for gas as 'dark' and infer that it amounts to an increase in 
the atomic gas by 28% (though due to opacity effects this could be as much as 
50%) and that the molecular hydrogen mass could double. On their own these 
values are not large enough to account for galactic rotation, but when combined 
with the uncertainties due to the low metalicty of gas in the outer regions of 
galaxies they become interesting. 

Although the emission from cold (T < 100 K) molecular hydrogen is rather 
weak there has been a detection of emission lines from it. Using ISO Valentijn 
and van der Werf (1999) detected H2 emission at 17 and 28^m along the plane 
of the edge on galaxy NGC891. This emission they assign to a warm (130-200K) 
and colder (80-90K) molecular component, with the colder component becoming 
more dominant at larger distance from the galaxy's centre. They conclude that 
the H2/HI ratio is 3 in the inner regions, but that this rises to greater than 20 
in the outer regions (fc = 21). Again this is very close to the value required to 
account for galactic rotation. 

So, if recent measures of molecular hydrogen indicate that there is very much 
more of it around than previously thought, what about the convincing three in- 
dependent measures of Xco described above? The first general comment is that 
the universally constant H2 CO relation is almost certainly far too simplistic 
possibly because all three measures were derived using things quite local to the 
Sun. Considering each measure in turn: 

1. There is an assumption here of a fixed gas-to-dust ratio throughout the 
galaxy, but it probably changes with the metalicity gradient (Magrini et 
al. 2011). 

2. The interpretation of the 7 ray data is very model dependent and critically 
dependent on both the local cosmic ray flux and gas density and how both 

^Drain ct al. (2007) use a value of Xco = 4.0 X 10^" cm-^ K'^ km"! s. 
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of these change throughout the galaxy. Pfennigcr et al. (1994) point out 
that the exponential scale length of 7 ray emission from the Galaxy is 
about three times longer than that of the sources that produce it (cosmic 
rays from pulsars and massive stars?). They suggest that the 7 rays are 
actually tracing '..some hidden form of hydrogen'. This '7-ray problem' 
has not gone away being recently confirmed by Ackermann et al. (2011). 
Strong ct al. (2004) suggest that the solution to this problem is that the 
value of Xco should be increased by a factor of from 5 to 10 when going 
from the inner to outer parts of the Galaxy. 

3. Molecular masses derived from the virial theorem are of course only appli- 
cable to gravitationally bound clouds. A problem for our cold gas 'idea' 
for the 'missing' mass is why does it not form stars? We do not have a 
definitive answer for this except that it maybe a threshold effect (Kenni- 
cutt 1989), but the threshold is just a bit higher than previously thought. 
If the gas is not forming stars then it is probably not in the form of gravi- 
tationally bound clouds and the logic of this method of obtaining H2 from 
CO no longer applies. 

4.3 Comments on the gas mass of galaxies 

From the above discussion we believe that it is entirely plausible that galaxies 
contain significantly more gas than has been previously thought. This plausible 
gas could easily have sufficient mass [k « 27) to account for the rotation of 
galaxies. Both cold atomic and molecular gas can be easily hidden in a galaxy 
and the observation of cold atomic hydrogen in our galaxy adds weight to our 
'missing' or 'dark' gas hypothesis (Knee and Brunt 2001). 

Some comments on the formation of one form of hydrogen from the other 
may put things in perspective. As stated earlier the THINGS team assume that 
the almost constant observed column density of HI across the discs of galaxies 
is due to the formation of molecular hydrogen. Molecular hydrogen is thought 
to form by surface recombination on dust grains. The process depends on the 
flux of atomic hydrogen, the number density and cross-section of the grains and 
a poorly known 'H2-formation efficiency factor'. Given the observed metalicity 
and corresponding gas-to-dust gradients (Magrini et al. 2011) across galaxies 
one might infer that molecular gas formation is suppressed in the outer regions 
of galaxies and hence here is the realm of the atom. 

H2 is also destroyed by photodisassociation in the interstellar medium unless 
it is self-shielded in dense clouds. Allen (2002) discuss an alternative view of 
these processes. He suggests that hydrogen in the discs of galaxies is primarily 
in the molecular form and that the atomic hydrogen we detect through the 21cm 
line is almost entirely the consequence of this photodisassociation. He goes on 
to model this process showing that the almost constant observed column density 
of atomic hydrogen is successfully explained in this way. This is a 'chicken and 
egg' situation - which came first the atomic or the molecular form? 

We now have three ways of explaining the almost constant surface density of 
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atomic hydrogen, it is cither; the result of the onset of molecule formation, an 
opacity effect and/or the result of the disassociation of molecules. These three 
possibilities reflect our uncertainty in the relation between atomic and molecular 
hydrogen in galaxies and hence the possibility that more of both exist. 

Finally, can big bang nucleosynthesis (BBNS) constraints on the baryon den- 
sity provide sufficient material to account for the rotation of galaxies? Pagel 
(1999) provides a nice summary of the baryon accounting and its relation to 
BBNS. From BBNS he constrains the density of baryons to be 0.035h^Q. Quot- 
ing Fukugita et al. (1998) Pagel gives a density in stars of O.OOSS/ifQ^ a factor 
of 10 lower than BBNS. When hot gas in clusters is included this raises the 
'detectable' baryon density to O-OOQlh^Q, now a factor of about 6 different from 
BBNS. This discrepancy between observed and predicted baryons is typically 
accounted for by invoking a huge reservoir of baryons in the warm intcr-galactic 
medium (see the introduction of Takei et al. 2011 for a recent discussion of the 
issue). An alternative is that they actually reside in the disc of galaxies. In 
section 3 we showed that the observational data requires a factor of about 7 
increase in the baryonic mass of galaxies to be consistent with their dynamics, 
which we can now see is in good agreement with the early universe constraints 
on the baryon mass density. This additional 'puzzle' should be added to our list 
in the introduction: 'The fractional increase in the mass of a galaxy required 
to account for its rotation is just that needed to account for the discrepancy 
between baryons observed and those inferred from big bang nucleosynthesis'. 

5 Other relevant issues 
5.1 An exponential disc of stcirs 

Just why the stars in a disc galaxy should have an exponential surface brightness 
distribution (Freeman, 1970) has been a problem for both those who assume 
disc galaxies form by monolithic gravitational collapse (Seiden et al. 1984) 
and those who more recently accept their hierarchical assembly (Button, 2009). 
In his seminal paper Freeman (1970) says that the answer to the exponential 
nature of discs is not known, but quoting Toomre (1964) he suggests that their 
origin may be something to do with the radial velocity dispersion of the stars, 
which needs to exceed a critical value to provide stability. Seiden et al. (1984) 
argue that a power law (Mestel?) disc is the most likely form of the gas disc 
and that an 'exponential like' surface density distribution for the stars comes 
aboTit from the star formation prescription they use - Stochastic Star Formation 
(SSF). Although SSF may offer a satisfactory explanation of the way stars form 
in irregular, particularly dwarf irregular, galaxies it is not so clear that it can 
explain star formation in galaxies like the Milky Way. It is quite clear that 
star formation in spiral galaxies is concentrated along the spiral arms and so 
the distribution of stars must logically be related to the natiirc of the spiral 
density wave. Button's (2009) paper highlights the fact that after at least 40 
years this issue is still not understood. The frame work for his paper is how to 
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get exponential discs from numerical simulations of the hierarchical assembly of 
galaxies, a disc formation mechanism very different to what Mestel envisaged 
and what is discussed here. 

It is rather straight forward to sec how the distribution of stars may be 
quite different to the distribution of the gas. For a Schmidt star formation 
law (Kennicutt, 2008) in which the Star Formation Rate per unit area (SFR) 
is proportional to some power of the gas surface density (SFR oc S^^^ where 
n ~ 1.4) the stars and gas will have power law surface density slopes of — n and 
-1 respectively. This simply assumes that the gas is distributed as in Mestel's 
disc and that star formation continues for the same length of time at all places 
in the disc. The real situation is obviously much more complex with the local 
stellar surface brightness dependent on the local physical conditions of the gas 
(propagation of a spiral density wave, for example) and the age of the stars. 

We conclude that the exponential nature of the surface brightness distri- 
bution of stars in disc galaxies provides no strong constraint on the overall 
distribution of the baryons because the star formation process is related to the 
gas surface density in a non-linear way. 

5.2 Measuring the local disc mass of our galeLxy 

If the majority of the mass of galaxies is confined to the disc and not the halo 
then it should be apparent in derivations of the disc surface density. One way of 
measuring the local surface density of the Milky Way disc is to measure the scale 
height and velocity dispersion of stars. Bahcall (1984) did this and found a value 
of about 83 Mq pc^, Kuijkcn and Gilmore (1991) 71 Mq pc^, and Holmberg 
and Flynn (2004) 56 Mq pc^. With a surface density of about 10^^ cm~^ at 
10 kpc from the galactic centre our 'typical' model galaxy is well within these 
values at 76 Mq pc^. If we use a simple exponential model for the distribution 

2 

of the stars above the plane then the predicted scale height is: /? = j^q^ where 
cr is the velocity dispersion of the stars. With cj = 30 km s~^ and S = 10^^ 
atoms cm~^ the predicted scale height is a reasonable 0.2 kpc. 

van der Kruit and Searle (1982) have also found that the scale height of 
the stars is reasonably constant with radius. They argue that this is due to 
heating of the stars by interactions with giant molecular clouds and the passage 
of the spiral density wave. If true the stars are not actually particularly good 
test particles for measuring the surface mass density of the disc. But another 
possibility for the constant scale height is that the disc surface mass density, as 
in Mestel's disc, falls much more slowly with radius than the light from the stars 
would predict (Fig. 10). Oiling (1995) has used the flaring of the HI disc at large 
galactic centre distances to show that the dark matter must reside in a flattened 
structure with axial ratio similar to a disc. In support of this Kalberla et al., 
(2007) say (they are modelling the atomic gas distribution of the Milky Way) 
that '..it is much easier to match a model to the observations if one considers 
dark matter associated with the Galactic disk.' 
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5.3 Disc stability 

That galaxies require massive halos, in addition to a disc, to give them stabihty 
has been another often quoted reason for invoking their existence. The origin of 
this idea can be traced back to a paper by Ostriker and Peebles (1973). They 
were carrying out a study of the stability of galactic discs using a numerical 
simulation. They found that models of the disc of our galaxy were 'grossly' 
unstable to bar-like modes on very short timescales, but that models with, in 
their case a rigid, massive halo, were much more stable. They argued that a 
halo with a mass from 1-2 times that of the disc would provide the necessary 
stability, this is much smaller than the size of the halo now required to produce 
the observed rotation. Ostriker and Peebles (1973) suggested a stability criteria, 
based on the ratio of kinetic and potential energy of the disc and halo, that could 
be used to decide if a stable disc forms. Criteria like this have since been used 
by numerical modellers of galaxy populations so that they can identify in their 
codes dark matter halos that are firstly bound and then those that will form 
different types of galaxies - for example bar unstable galaxies are surprisingly 
often labelled as elliptical. It is not clear how useful these global criteria are 
in deciding whether a gas disc is actually stable or not within a simulation, or 
whether this is relevant to how a real galaxy evolves with time (Athanassoula, 
2008). For example Ostriker and Peebles (1973) suggested a massive halo as a 
means to suppress bar formation yet we now know that a large fraction of disc 
galaxies actually do have bars of one form or another. It is also quite clear from 
observations of the gas discs of galaxies that instability plays an important role 
in the way galaxies evolve and particularly how they form stars. Many disc 
galaxies are observed to have, along with bars, a large number of other distinct 
structiires that form in the gas - spiral arms, filaments and clumps. If galaxies 
do have massive dark halos then they are not actually doing a very good job of 
preventing instabilities in the disc. 

Ostriker and Peebles (1973) actually showed that the discs they modelled 
became stable when the ratio of kinetic energy of rotation to gravitational energy 
was reduced to f « 0.14. As this ratio for a 'typical' galaxy was higher than 
this they favoured a halo as a solution. For our 'typical' galaxy t w 0.5, so at 
face value this disc is unstable because it is not hot enough. 

Efstathiou et al. (1982) questioned the validity of the Ostriker and Peebles 
stability criteria and suggested a new test: 



where Vm is the maximum velocity of the disc, Md is the disc mass and Rd 
the exponential disc scale length. It is difficult to relate this to the Mestel disc 
as the criteria was developed for an exponential disc. However, t* looks like 
the square root of the kinetic energy per unit mass divided by the square root 
of the potential energy per unit mass in which case the value of t* would be 
about 1.4 for the our 'typical' galaxy and the stability criteria would actually 
be satisfied. Subsequently Athanassoula (2008) has shown that this criterion 
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{t* > 1.1) is also flawed and that it docs not consistently in simulations, predict 
which galaxies will or will not produce strong bars. 

In summary it is not clear if any of these stability criteria tell us very much 
about whether a disc will maintain its gross propcritics over many rotation 
periods, which is our basic requirement. So, we do not believe that the stability 
of the disc provides a strong constraint on the necessity for an extensive massive 
halo (see also p. 603 Binney and Tremaine, 1994). 

The Toomre instability (Toomre, 1964) is most often quoted as the means 
by which stars form, the criteria being: 



where k is the cpicyclic frequency and c is the velocity dispersion of the gas 
(taken to be the sound speed c = \ r^)- For a Mestel disc: 



which interestingly is constant with radius and equal to 0.02 for our 'typical' 
disc. So the prediction is that the disc will fragment and form stars, just as we 
would want it to do. Toomre stability requires a much hotter disc of about 10^ 



If discs are to form in the way Mestel envisaged then initially there must 

have been an adiabatic collapse; which gave rise to the gross features of the gas 
disc. This seems a reasonable assumption because in the early Universe we 
expect the gas to be of low metalicity and so there will be few metals to provide 
the necessary cooling lines. Subsequent radiative processes in metal enriched 
gas can produce the smaller scale structures and allow the gas to cool and form 
stars as the galaxy evolves. This is just the scenario proposed by those who 
first considered the collapse and siibsequcnt fragmentation of proto-galactic gas 
clouds (Hoyle, 1953, Silk, 1977). The critical parameter (Silk, 1977) is the ratio 
of the cooling time [tcooi oc where p is the density) to the gravitaional 
collapse time {torav oc p~^^'^). For tcooi/tcrav oc p~^^'^ >> 1 the collapse is 
adiabatic, but as the collapse continues, the density increases and the cooling 
time becomes shorter, eventually because of the nature of the hydrogen gas 
cooling curve the gas becomes isothermal and begins to fragment into clouds 
that eventually form stars. 

5.4 The gas consumption problem 

The average timescale for gas consumption by star formaton in our galaxy is 
of order 4 Gyr (Larson et al., 1980). From observations of the age, metalicity 
and velocities of local stars it appears that the star formation rate in the disc 
of our galaxy has been almost constant for billions of years longer than this. 
This almost constant star formation rate has recently been confirmed not just 
for our galaxy, but also for the Local Group as a whole by a detailed study of 





K. 
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stellar populations (Drozdovsky ct al., 2008). Surely as the gas is depleted the 
star formation rate should decrease? Larson et al.'s explanation, and that now 
normally accepted, is that there is continual infall of fresh gas into the galactic 
disc. The problem is that significant amounts of in falling gas with varying 
velocity and direction will thicken the galactic disc beyond what is observed 
(Quinn, 1987, Toth and Ostriker, 1992). This does not need to be the case if 
the stars already sit in a huge reservoir of gas in which they only represent a 
small fraction of the mass. 

5.5 Constraints imposed by quasar absorption lines 

With typical atomic hydrogen column densities of lO^^^^'^ cm^^ (Fig. 10) sight 
lines to quasars through our model 'typical' disc would produce damped Ly^ 
absorption spectra. Observations of damped Ly^ absorption systems produce 
column densities at most 10^^ cm~^. At first sight this seems a problem for our 
massive Mestel disc hypothesis. However, from the reddening of background 
quasars Fall et al. (1989) have shown that there is dust associated with Lya ab- 
sorption features. Pettini et al., (1994) have shown that typical Lyo, absorption 
systems have metalicitics about 0.1 of the solar value. If the dust is depleted in 
the same way and has the same properties as dust in the Milky Way then we 
might expect Ay w ry = where this is just a scaling of the Burstein and 
Heiles (1978) relation in which they measure approximately a V band optical 
depth of unity through an HI column density {N^j) of 10^^ cm^^. Thus optical 
sight lines to quasars through gas at column densities greater than about 10^^ 
cm~^ will be optically thick and these quasars will be missing from optically 
selected catalogues (Wright, 1990, FaU and Pei, 1993). Webster et al., (1995) 
have argued exactly this point. They show that radio selected quasars are red 
at optical wavelengths compared to optically selected quasars. They put this 
down to optical extinction by dust and infer that there is a large population of 
radio quiet quasars that remain undetected because they are extinguished. The 
high frequency of damped Lya absorption features has normally been explained 
by assuming large extended HI halos at early epochs (Briggs and Wolfe 1983). 
This conclusion would need to be revisited if the actual column density through 
galactic discs is much higher than expected from the observed 21cm column 
densities of local galaxies. 

5.6 Dwarf galetxies 

A galaxy mass of M w lO'^ Mq and a size of _Ro ~ 5 kpc leads to a Mcistel disc 
surface density of Sq ~ 10^^ cm~^. Low mass galaxies like this must contain 
gas that is severely affected by the ionising background and no atomic hydrogen 
will be detected at 21cm. If the hydrogen is ionised then re-combination emis- 
sion lines, particularly E.^, should be observable in dwarf galaxies as a test of 
our hypothesis. This is a very similar experiment to that described by Bland- 
Hawthorn et al. (1997) in their search for ionised gas at the edge of nearby 
galaxies. 
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Wc know that at least back to a rcdshift of about z = 1 the star formation 
rate per unit comoving volume was higher than it is today (Madau et al. 1998). 
As the Universe was also smaller the energy density of the ionising radiation 
must have been considerably higher. So, as the Universe ages from z » 1 the 
ionising flux decreases and the lower column densities of atomic hydrogen in low 
mass galaxies can resist ionisation. The HI can then cool and eventually form 
stars. If lower mass galaxies have lower column density discs (E cx Ad^^"^) then 
there should be a relationship between their redshift and the ages of their stellar 
populations (downsizing?). We may also speculate that early on in the history 
of the Universe, before the majority of stars were formed, the ionising flux would 
again be low and all galaxies would be able to form stars, thus accounting for 
the old stellar population found in all dwarfs. 

It is observed that brighter dwarf galaxies have the lowest mass-to-light 
ratios something that is explainable by the baryons they contain (Mateo 1998). 
It is only for the lowest luminosity dwarf galaxies that the mass-to-light ratios 
become extreme, explained here by the onset of ionisation. Thus we suggest that 
there are three regimes: big (massive) galaxies appear to have high mass-to-light 
ratios because they have an as yet undetected cold gas component, intermediate 
galaxies have mass-to-light ratios fully accounted for by their baryons, small 
galaxies have high mass-to-light ratios because they are ionised. The boundaries 
being roughly at 10^ and 10^ Mq. 

6 The rotation curve puzzles and TuUy-Fisher 
revisited 

With the optically thick Mestel disc illustrated in Fig 10 it is quite easy to see 
that the observed HI surface density is very close to a scaled version of the real 
surface density of hydrogen, thus explaining Bosma's (1978, 1981) observation 
that the rotation predicted from the gas is just a scaled version of the total. We 
are not suggesting that the gas is totally opaque and atomic, but that there could 
plausibly be very much higher surface densities of both atomic and molecular 
hydrogen than have previously been thought. 

For our hypothesised massive baryonic disc the largest fraction of the mass 
of a galaxy is now in the gas so there is no longer any reason to use maximum 
disc models based on the distribution of the stars. There is also now no longer 
any disc halo conspiracy because there is no longer any need for a halo. 

In the introduction we used a dynamical argument to try to explain the 
TuUy-Fisher relation, this predicted a slope of 1/4. In fact the observed slope 
is closer to 1/3. Binney and Tremaine (1987) give exactly a value of 1/3 for 
the B band Tully-Fisher relation in precise agreement with what we found in 
section 3 (Fig. 7). Mestel shows that to form his disc via collapse of a spherical 
gas cloud, the density of the cloud (p) must be approximately constant. If 
this is true and about the same density for all galaxies then all Mestel discs 

1 /3 

will satisfy the Tully-Fisher relation Vm oc Mj!^^ (with luminosity as a good 
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proxy for baryonic mass). Mcstcl discs are predicted to follow the Tully-Fisher 
relation with logarithmic slope of 1/3 because of the structure of the proto- 
galactic clouds they formed from. No conspiracy between surface brightness 
and mass-to-light ratio is required. 

7 Summary 

The observed rotation curves of galaxies do not necessarily imply that galaxies 
contain non-baryonic dark matter residing in an extended halo. The collapse of 
baryonic constant density/isothermal spheres leads naturally to surface density 
distributions that explain galactic rotation as long as one assumes that there is 
about a factor of 7 more baryonic material in the disc than is currently accounted 
for. The observed almost constant surface density of atomic hydrogen may be 
explained by a combination of cold opaque atomic gas with additional cold 
molecular hydrogen and may be sufficient to account for this missing factor of 
7. 

The massive baryonic disc has many attractive features. All 'conspiracies' 
between the spherical dark matter halo and the baryonic disc are removed, it 
provides a reservoir of gas for continued star formation over a Hubble time 
(something possibly already identified in the 7-ray data) and it removes the 
necessity for unknown dark matter. However, currently there is no conclusive 
evidence for so many baryons in the disc and future observations of atomic and 
molecular hydrogen in galaxies are required to justify or refute our hypothesis. 
This though should be contrasted with the invention of dark matter and the 
continued lack of direct evidence for its existence. It should also be a concern 
that the additional baryonic mass required to account for rotation is just that 
required to reconcile observed baryons with those predicted by BBNS. 

This paper does not try in any way to address other observations that im- 
ply the existence of dark matter, such as baryon acoustic oscillations, peaks 
in the cosmic microwave power spectrum and studies of large scale structure, 
though there are some who suggest that these too can be explained in other 
ways (Shanks, 2005). Here we only address the issue of the interpretation and 
implications of galaxy rotation curves. However we note in passing that other 
problems with the current ACDM model may have more straight forward solu- 
tions as a result: 

1. Not enough low luminosity galaxies - if the major formation process is 
by monolithic collapse then this problem goes away as the mass spectrum 
of galaxies is a faithful reflection of the density fluctuations that created 
them. 

2. The cusp core problem - no dark matter halo so no problem. 

3. Too many bright galaxies predicted by ACDM - see (i). above. 

4. The 'bias' problem is removed as baryons trace the mass. 
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Wc conclude that it may be more productive to look for a dark baryonic com- 
ponent associated with hydrogen than to invoke and search for candidate exotic 
dark particles that have no other reason to exist other than to satisfy a whim 
or WIMP. 
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